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ABSTRACT 

We present the Suzaku broad band observations of two AGNs detected by the Swift/BAJi hard X-ray (>15 
keV) survey that did not have previous X-ray data, Swift J0601.9-8636 and Swift J0138.6-4001. The 
Suzaku spectra reveals in both objects a heavily absorbed power law component with a column density of 
A^h — 10 23 5 ~ 24 cm" 2 that dominates above 10 keV, and an intense reflection component with a solid angle £ 2tt 
from a cold, optically thick medium. We find that these AGNs have an extremely small fraction of scattered 
light from the nucleus, 0.5% with respect to the intrinsic power law component. This indicates that they 
are buried in a very geometrically-thick torus with a small opening angle and/or have unusually small amount 
of gas responsible for scattering. In the former case, the geometry of Swift J0601. 9-8636 should be nearly 
face-on as inferred from the small absorption for the reflection component. The discovery of two such objects 
in this small sample implies that there must be a significant number of yet unrecognized, very Compton thick 
AGNs viewed at larger inclination angles in the local universe, which are difficult to detect even in the currently 
most sensitive optical or hard X-ray surveys. 

Subject headings: galaxies: active — gamma rays: observations — X-rays: galaxies — X-rays: general 



1, INTRODUCTION 

Many observations imply the presence of a large number of 
heavily obscu red Active Galactic Nuclei (AGNs) in th e local 
universe (e.g.. lMaiolino et al.|[l99llRisaliti et al.lll999l) . The 
number density of AGNs subject to absorption with a line- 
of-sight Hydrogen column density log A^h ^ 23.5 cm" 2 is a 
key parameter in understanding the accretion his tory of the 
universe (e.g.. lSetfi & Wo ltier 19891: lFabian~lll999l) . Accord- 
ing to synthesis models of the X-r ay background (XRB) (e.g., 
lUeda et aLll2.003b iGilli et al.l2007l) . much of the peak intensity 
of the XRB at 30 keV should be produced by these objects. 
In spite of the potential importance of their contribution t o 
the growth of supermassive black holes (Mar coni et al.l2004l) . 
the nature of this population of AGNs, even in the local uni- 
verse, is only poorly understood due to strong biases against 
detecting them. In these objects, the direct emission in the 
UV, optical, and near IR bands as well as at E < 10 keV from 
the nucleus is almost completely blocked by obscuring matter, 
making it difficult to probe the central engine. 

Hard X-ray surveys at energies above 10-15 keV provide us 
with an ideal opportunity to select this population of AGNs 
as long as the column density is less than log A^h ~ 24.5 
cm" 2 . Recent surveys perf ormed with Swi/f/BAT (15-200 
keV:lMarkwardtetal 1120051) and INTEGRAL (10-100 keV; 
Bassa ni et alj 12006: Beckmann et al. 2006; ISazonov et al.l 
2007), because of their relative insensitivity to absorption, are 
providing one of the most unbiased AGN samples in the lo- 
cal universe including Compton thick AGNs, i.e., those with 
log A^h > 24 cm" 2 . In fact, these surveys have started to de- 
tect hidden AGNs in the local universe located in galaxies that 
were previously unrecognized to contain an active nucleus at 
other wavelengths. 
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To unveil the nature of these new hard X-ray sources, 
follow-up observations covering a broad energy band are cru- 
cial. In this paper, we present the first result s from follow- 
up observations with the Suzaku observatory ( Mitsu da et al.1 
[2007b of the AGNs Swift J0601.9-8636 and Swift J0138.6- 
4001 detected in the SwiftlBPS survey. These targets are es- 
sentially randomly selected from a bright Swift AGN sample 
for which soft X-ray (<10 keV) spectroscopic observations 
had never been performed, and are thus reasonable represen- 
tatives of unknown AGN populations selected by hard X-rays 
above 15 keV. In § 3, we also present an optical spectrum 
of Swift J0601. 9-8636 taken at the South African Astronom- 
ical Observatory (SAAO) 1.9-m telescope. The cosmologi- 
cal parameters (H tt m , tt x ) = (70 km s" 1 Mpc" 1 , 0.3, 0.7) 
(Spergel et al. 2003) are adopted throughout the paper. 

2. THE SUZAKU OBSERVATIONS AND RESULTS 

Table 1 summarizes our targets and observation log. Swift 
J060 1. 9-8636 is opti cally identified as the galaxy ESO 005- 
G004 (lLaubertsl["l982l) at z = 0.0062 with no previous firm 
evidence for AGN activity. The optical coun terpart of Swift 
J0138.6-W01 is the galaxy ESO 297-G018 (lLaubertslll982h 
at z = 0.0252, which was identified as a narrow line AGN 
dKirhakos & Steinerll 19901) . 

Suzaku, the 5th Japanese X-ray satellite, carries four sets 
of X-ray mirrors each with a focal plane X-ray CCD camera, 
the X-r ay Imaging Spectrom eter (XIS-0, XIS-1, XIS-2, and 
XIS-3: iKovama et all 120071) . and a n on-imaging i nstrument 
called the Hard X-ray Detector (HXD: lTakahashi et al.ll2007l) . 
which consists of the Si PIN photo-diodes and GSO scintil- 
lation counters. The XIS and PIN simultaneously covers the 
energy band of 0.2-12 keV and 10-70 keV, respectively. The 
unique capabilities of Suzaku, high sensitivity in the 12-70 
keV band and broad band coverage with good spectral reso- 
lution, are critical for studies of highly absorbed AGNs. 

We observed Swift J0601.9-8636 and Swift J0138.6-4001 
with Suzaku on 2006 April 14 and June 5 for a net exposure 
of 20 and 21 ksec (for the XIS), respectively. Standard anal- 
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TABLE 1 
Targets and Observation Log 



Swift 


Optical ID 


redshift 


Start Time (UT) 


End Time 


Exposure' 1 


J060 1.9-8636 
J0138.6-4001 


ESO 005-G004 
ESO 297-G018 


0.0062 
0.0252 


2006/04/13 16:24 
2006/06/04 18:13 


04/14 01:52 
06/05 05:00 


19.8 ksec 
21.2 ksec 



a Based on a good time interval for the XIS-0. 



ysis was made on data products, which were processed with 
the latest calibration (version 1.2). We detected both sources 
at high significance with the XIS and PIN. The XIS spectra 
were accumulated within a radius of 2 arcmin around the de- 
tected position. The background was taken from a source free 
region in the field of view. The spectra of three Front-side Illu- 
minated CCD (hereafter FI-XIS; XIS-0, 2, and 3) are summed 
together, while that of Back-side Illuminated CCD (BI-XIS; 
XIS-1) is treated separately in the spectral fit. Examining the 
spectra of 55 Fe calibration source, we verify that the energy 
scale and resolution is accurate better than 10 eV and 60 eV 
levels, respectively, at 5.9 keV. For the analysis of the PIN, 
we utilized only data of Well units with the bias voltage set 
at 500 eV (W0,l,2,3 for Swift J0601.9-8636 and Wl,2,3 for 
Swift J0138.6-4001) with the best available models of PIN 
background provided by the HXD team 5 . 

To obtain the best constraint from the entire data, we per- 
form a simultaneous fit to the spectra of XIS (FI and BI), 
PIN, and the archival Swift BAT, which covers the 0.2-200 
keV band as a whole. The BAT spectra consist of four en- 
ergy bins over the 15-200 keV range and are useful to con- 
strain the power law index. Here we allow the relative flux 
normalization between Suzaku and Swift (BAT) to be a free 
parameter, considering possible time variability between the 
observations. We fixed the normalization ratio between the 
FI-XIS and the PIN based on the calibration result using the 
Crab Nebula. All the absolute fluxes quoted in this paper refer 
to the flux calibration of the FI-XIS. We find that the 10-50 
keV PIN fluxes of Swift J0601.9-8636 and Swift J0138.6- 
4001 are 1 x 10" 11 and 4 x 10~ n erg cm" 2 s" 1 , indicating time 
variability by a factor of 0.5 and 1.6, respectively, compared 
with the averaged flux measured by the Swift/BAT over the 
past 9 months (Tueller et al, in preparation). 

Figure 1 shows the FI-XIS and PIN spectra unfolded for the 
detector response (for clarity the BI-XIS and BAT spectra are 
not plotted). The X-ray spectrum of Swift J0601. 9-8636 be- 
low 10 keV is dominated by a hard continuum with few pho- 
tons below 2 keV, consistent with the previous non-detection 
in soft X-rays (an upper limit of 1 .3 x 10~ 13 erg c m" 2 s" 1 in the 
0.1-2 .4 keV band by the RO SAT All Sky Survey: I Voges etalj 
2000). We find that the broad band spectrum can be well re- 
produced with a model consisting of a heavily absorbed power 
law with log A^h — 24 cm" 2 , which dominates above 10 keV, 
and a mildly absorbed reflection component from cold mat- 
ter accompanied by a narrow fluorescence iron-K line, which 
dominates below 10 keV. The large column density is consis- 
tent with the ob served equivale nt width (EW) of the iron-K 
line, w 1 keV dLevenson et alJl2002h . Swift J0138.6^001 
shows a similar spectrum but with a smaller absorption of 
log A^h = 23.7 cm" 2 for both transmitted and reflected com- 
ponents. 

5 verl.2_d for Swift J0601.9-8636 and verl.2_wl23 for Swift J0138.6- 
4001 
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FIG . 1 .— The broad-band energy spectra of (a) Swift J0601 .9-8636 and (b) 
Swift J0138. 6^1001 unfolded for the detector response in units of E 2 F(E), 
where F(E) is the photon spectrum. For clarity, we only plot the summed 
spectrum of the three FI-XIS (below 12 keV), and those of PIN (above 12 
keV), while the spectral fit is performed to the whole XIS+PIN+BAT data 
including the BI-XIS. The crosses (black) represent the data with Icr statisti- 
cal errors. The histograms show the best-fit model with separate components. 
The upper solid line (red), dot-dashed line (light blue), lower solid line (blue), 
and dotted line (magenta) correspond to the total, iron-K emission line, re- 
flection component, and scattered component, respectively. 

The spectral model is represented as 

F(E) = e-° {E)N ™ [fAE- r +e-° (E)N »AE- r +e-° (E)N '» a C(E)+G(E)] , 

where A^ al is the Galactic absorption column density fixed 
at 2.0 x 10 20 cm" 2 for both targets dDickev & Lockmanl 19901) . 
A^h the local absorption column density at the source redshift 
for the transmitted component, that for the reflected com- 
ponent (assumed to be the same as A^h for Swift J0138.6- 
4001), and cr(E) the cross section of photo-electric absorption. 
The term C(E) re presents the reflection compo nent, calculated 
using the code in Magdziarz & ZdziarskJ( ll995h : we leave the 
solid angle 57 of the reflector as a free parameter by fixing the 
inclination angle at 60° and cutoff energy at 300 keV, assum- 
ing Solar abundances for all elements. R(= $7/271") > 1 means 
the transmission efficiency should be ^ l/R (see § 4). The 
term G(E) is the narrow iron-K emission line modelled by 
a Gaussian profile, where we fix the lcr width at 50 eV to 
take account of the response uncertainty (and hence, the line 
should be considered to be unresolved). The best-fit param- 
eters are summarized in Table 2 with the intrinsic 2-10 keV 
luminosity, L2-10, corrected for the absorption and transmis- 
sion efficiency of l/R. 

We confirm that these results are robust, within the statis- 
tical error, given the systemati c errors in the back ground es- 
timation of the PIN detector dKokubun et al.1 120071) . In the 
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TABLE 2 
Best Fit Spectral Parameters 





Swift 


J0601. 9-8636 


J0138.6^1001 


(1) 


N H (10 22 cnr 2 ) 


101 1s 


46±4 


17) 
(3) 
(4) 


p 
R 

Nfi a (10 22 cnr 2 ) 


1 qc+0.36 
1 7+ 3 - 5 


loD -0.04 
7 1+0.4 
z - 1 -1.2 

<=m 


(5) 


/scat (%) 


0.20 ±0.11 


23+ 023 


(6) 


£cen (keV) 


6.38 ±0.02 


6.38 ±0.03 


(7) 


E.W. (keV) 


1.06 ±0.16 


0.20 ±0.05 


(8) 


F 2 -10 (ergcnT 2 s -1 ) 


1.1 x 10~ 12 


3.3 x 10~ 12 


(9) 


F 10-50 (erg cm" 2 s" 1 ) 


9.8 x 10~ 12 


3.9 x 10-" 


(10) 


L2-10 (erg s" 1 ) 


8.3 x 10 41 


3.9 x 10 43 




X 2 (dof) 


20.0 (27) 


101.3 (87) 



NOTE. — (1) The line-of-sight hydrogen column density for the transmit- 
ted component; (2) The power law photon index; (3) The relative strength of 
the reflection component to the transmitted one, defined as R = Q/2tt, where 
Q is the solid angle of the reflector viewed from the nucleus; (4) The line- 
of-sight hydrogen column density for the refection component (assumed to 
be the same as Nh for Swift J0138. 6^1001); (5) The fraction of a scattered 
component relative to the intrinsic power law corrected for the transmission 
efficiency of l/R when R > 1; (6) The center energy of an iron-K emission 
line at rest frame. The lcr line width is fixed at 50 eV; (7) The observed equiv- 
alent width of the iron-K line with respect to the whole continuum; (8)(9) 
Observed fluxes in the 2-10 keV and 10-50 keV band; (10) The 2-10 keV 
intrinsic luminosity corrected for the absorption and transmission efficiency 
of l/R. 

The errors are 90% confidence limits for a single parameter. 

case of Swift J0138. 6^4001, where photon statistics is domi- 
nated by the XIS data, we limit the allowed range of photon 
index to T = 1.63-2.02 in the simultaneous fit, being con- 
strained from the BAT spectrum. We have limited R < 2.5 
for this source to avoid physical inconsistency between R and 
the EW of an iron-K line; otherwise (i.e., in more "reflection- 
dominated" spectra), we should expect a larger EW than the 
observed value of « 0.2 keV. 

3. OPTICAL SPECTRUM OF SWIFT J0601.9-8636 

We performed an optical spectroscopic observation of 
Swift J060 1.9-8636 (ESO 005-G004) during the night of 
2007 March 16, using the SAAO 1.9-m telescope with the 
Cassegrain spectrograph. Grating six, with a spectral range 
of about 3500-5300 A at a resolution around 4 A, was used 
with a 2 arcsec slit placed on the center of the galaxy for 
a total integration time of 2400 s. To derive the sensitivity 
curve, we fit the observed spectral energy distribution of stan- 
dard stars with low-order polynomial. The co-added, flux- 
calibrated spectrum in the 4000-5500 A range is shown in 
Figure 2. It reveals a rather featureless spectrum with no ev- 
idence for H/3 or [O III] A5007 emission lines, typical for 
this type of non-active edge-on galaxy within this spectral 
range. The 90% upper limit on the [O III] flux is conserva- 
tively estimated to be 3 x 10~ 15 erg cm" 2 s" 1 , corresponding 
to a luminosity of 3 x 10 38 erg s" 1 . This yields the ratio of 
the intrinsic 2-10 keV luminosity to the observed [O III] lu- 
minosity of > 2800. Although the (unknown) extinction cor- 
rection for [O III] could reduce the value, the result is con- 
sistent with Swift J0601. 9-8636 having an intrinsically weak 
[O III] emission relative to hard X-rays compared with other 
Seyfert galaxies dBassani et al.|[l999l: iHeckmann et al.l 120051; 
Netz er et alJl2007ir In particular, this object would not have 
been selected to be an AGN on the basis of its [O III] or H/3 
emission. 

4. DISCUSSION 
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FIG. 2. — The optical spectrum of the nucleus region of Swift 10601.9- 
8636 (ESO 005-G004) in the 4000-5500 A wavelength range, taken with the 
SAAO 1.9-m telescope. The arrow denotes the position of the [O III] A5007 
line. 

Both sources show an intense reflection component relative 
to the transmitted one. Using the standard reflection model 
( Magdziarz & Zdziarski 1995), we find that the solid angle 
of the reflector n/2ir viewed from the nucleus exceeds unity, 
which is apparently unphysical if attributed only to geome- 
try. This implies that a part of the direct emission is com- 
pletely blocked by non-uniform material in the line of sight 
even above 10 keV. The reflection-dominated nature of the 
spectra of heavily obscured AGNs, if common, has an impact 
on the population synthesis model of the XRB, wh ere a much 
weak er reflection is assumed for type 2 AGNs dGilli et alj 
120071) . Another possibility is that this apparent very high re- 
flection fraction is due to time variability, that is, the decrease 
of the flux in the transmitted light is echoed with a time de- 
lay corresponding to the difference in light paths between the 
emitter, reflector, and observers. 

It is remarkable that both Swift J0601. 9-8636 and Swift 
J0138. 6-4001 have a very small amount of soft X-ray scat- 
tered emission, less than 0.46% of the intrinsic power law 
component. (If we fix = in the spectral fit of Swift 
J0601. 9-8636, then we obtain a photon index of 1.80 ±0.29 
and no significant scattered component with a 90% upper 
limit of 0.47%.) As far as we know these are amongst the 
lowest scattered fractions ever seen f rom an absorbed AGN 
dTurner et al.lll997t ICappi et al.l 120061) . In optically selected 
Seyfert 2 galaxies, the presence of prominent soft X-ray emis- 
sion is common (e.g., Guain azzi et al.|[2005b . Such emission 
probably originates from the same extended gas resp onsible 
for the optical [O III] emission dBianchi et al.l l2006h . This 
type of emission from "classical" Seyfert 2 galaxies has al- 
ways been seen in the spectra of objects well studied so 
far. However, this sample is dominated by optically selected 
Seyfert 2 galaxies, which require a scattered component to be 
selected. 

The scattered fraction is proportional to both the solid angle 
of the scattering region as viewed from the nucleus, f2 sca t, and 
the scattering optical depth, T scat . Hence, the observed small 
scattered fraction means small £7 scat and/or small r scat i.e., de- 
ficiency of gas in the circumnuclear environment, for some 
unknown reason. 

The first possibility, which we favor as a more plausi- 
ble case, indicates that these AGNs are buried in a very 
geometrically-thick obscuring torus. Assuming that the typi- 
cal scattering fraction of 3% corresponds to the effective torus 
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half-opening angle (see lLevenson et aT]|2002l for definition) 8 
of 45 degree, our results (< 0.5%) indicate 8 ^ 20 degree. In 
the case of Swift J0601. 9-8636, the small absorption for the 
reflection component, which probably comes from the inner 
wall of the torus, suggests that we are seeing this source in 
a rather face-on geometr y. Indeed, applying the formalism 
of iLevenson et al.l d2002l) to the observed EW of the iron-K 
line, we infer that the inclination angle with respect to the 
axis of the disk, i, is smaller than 40 degrees if 8 < 20 de- 
gree. For Swift J0138. 6-4001, the presence of a high column 
density for the reflection component implies a more edge-on 
geometry than in Swift J0601. 9-8636. The observed EW of 
0.20 ± 0.05 keV can be explained if the torus is patchy or 
has a geometrical structure such that the line-of-sight column 
density is much smaller than that in the disk plane. 

We infer that this type of buried AGNs is a significant frac- 
tion of the whole AGN population, although an accurate es- 
timate of this fraction is difficult at present due to the small 
number statistics 6 . The observed fraction of heavily obscured 
AGNs with log N n > 23.5 is ab out 25% among the ha rd X- 
ray (E > 15 keV) selected AGNs dMarkwardt et alj2005l) . The 
true number density of obscured AGNs could be much larger, 
however. If we saw the same system of our targets at much 
larger inclination angles (z ^> 40 degree), the observed flux 
of the transmitted component would be much fainter even in 
hard X-rays due to the effect s of repeated Compton scatter- 
ings (IWilman & FabiarJll999l) . Our results imply that there 
must be a large number of yet unrecognized, Compton thick 
AGNs in the local universe, which are likely to be missed even 
in the Swift and INTEGRAL surveys. 

The existence of AGNs wi t h a g eometrically thick torus 
was predicted by Fa bian et"ai1 d!998l) . where the extreme ob- 
scuration was postulated to be caused by a nuclear starburst. 
Using the 60/im and 100/im fluxes measured by Infrared As- 
tronomical Satelli te (IRAS), we obta in the far infrared lumi- 
nosity (defined by lDavid etaill 19921) of L F ir = 4.4 x 10 43 erg 
s" 1 and 7.1 x 10 43 erg s -1 , and hence the ratio between the 
2-10 keV to far infrared luminosities of La-io/Lpm ~ 0.02 
and re 0.5 for Swift J0601.9-8636 and Swift J0138.6-4001, 
respectively. While the result of Swift 10138.6^-001 is con- 

6 W e note that a similar object has also been found by Coma stri et alj 
120071) with the Suzaku follow-up of hard X-ray (> 10 keV) selected AGNs. 



sistent wit h those of the 2-10 ke V selected AGNs in the local 
universe dPiccinotti et al.llT982h within the scatter, the small 
L2-io/^fir ratio of Swift J0601. 9-8636 indicates a possibly 
significant starburst activity. However, this is not supported 
by the optical spectrum of this object, which apparently shows 
no evidence for a significant amount of star formation. The 
reason behind the difference between the two sources is un- 
clear. 

By using the unique combination of the Swift BAT survey 
and the Suzaku broad band spectral capabilities, we are dis- 
covering a new type of AGN with an extremely small scat- 
tering fraction. This class of object is most likely to con- 
tain a buried AGN in a very geometrically-thick torus. This 
population was missed in previous surveys, demonstrating the 
power of hard X-ray ( > 1 keV) surveys to advance our global 
understanding of the whole AGN population. In particular, 
we predict that the objects should have fainter [O III] emis- 
sion luminosity relative to the hard X-ray luminosity com- 
pared with classical Seyfert 2 galaxies because much less of 
the nuclear flux "leaks" out to ionize the narrow line gas. As 
shown above, the optical spectrum of Swift J0601. 9-8636 is 
consistent with this prediction. This study is particularly im- 
portant since the existence of numerous such objects would 
make surveys that rely on the [O III] emission incomplete by 
missing many of buried AGNs and incorrectly estimating the 
true AGN luminosity. 
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